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We investigate the behavior of flows, including turbulent flows, driven by a horizontal 
body-force and subject to a vertical magnetic field, with the following question in mind: 
for very strong applied magnetic field, is the flow mostly two-dimensional, with remaining 
weak three-dimensional fluctuations, or does it become exactly 2D, with no dependence 
along the vertical? 

We first focus on the quasi-static approximation, i.e. the asymptotic limit of vanishing 
magnetic Reynolds number Rm <C 1: we prove that the flow becomes exactly 2D asymp¬ 
totically in time, regardless of the initial condition and provided the interaction parameter 
N is larger than a threshold value. We call this property absolute two-dimensionalization: 
the attractor of the system is necessarily a (possibly turbulent) 2D flow. 

We then consider the full-magnetohydrodynamic equations and we prove that, for low 
enough Rm and large enough N, the flow becomes exactly two-dimensional in the long¬ 
time limit provided the initial vertically-dependent perturbations are infinitesimal. We 
call this phenomenon linear two-dimensionalization: the (possibly turbulent) 2D flow is an 
attractor of the dynamics, but it is not necessarily the only attractor of the system. Some 
3D attractors may also exist and be attained for strong enough initial 3D perturbations. 

These results shed some light on the existence of a dissipation anomaly for magneto¬ 
hydrodynamic flows subject to a strong external magnetic field. 


1. Introduction 

Liq uid m etal turbulence is encountered in various situ ations, rangin: 


( Davidson 1999) to the flow in the Earth’s outer core ( Moffatt 197 


from metallurgy 
, including labo¬ 
ratory experiments on magnetohydrodynamic (MHD) turbulence (lAlemanv et alJll97flt 
Gallet et al. HI). Because their kinematic viscosity is much lower than their mag¬ 
netic diffusivity, with magnetic Prandtl numbers Pm typically in the range of 10~ 6 to 
10 —5 , flows of liquid metals can be turbulent and still be characterized by a small mag¬ 
netic Reynolds number: in most laboratory experiments on MHD turbulence, the kinetic 
Reynolds number is in the range 10 4 — 10 6 , with a magnetic Reynolds number rarely 
exceeding unity. At small enough scales, turbulence in Earth’s outer core may also be 
considered as a low- Rm flow. 

Turbulence subject to a strong external magnetic field tends to become two-dimensional, 
the velocity field becoming weakly dependent on the coordinate along the applied mag¬ 
netic field (hereafter denoted as the vertical axis by convention). Indeed, any flow that 
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depends on the vertical coordinate shears the applied field and induces electrical cur¬ 
rents. For low magnetic Reynolds number, such currents are rapidly damped through 
Ohmic dissipation and the flow structures that depend on the vertical are more damped 
than vertically-invariant eddies. This intuitive argument can be made more precise in the 
asymptotic limit of vanishing magnetic Reynolds number Rm, the so-called quasi-static 
approximation. In this limit the applied fi eld simply appears as an additional Ohmic 
damp ing term in the Navier-Stokes equation ( Sommeria fe Moreau lll982 : Knaepen &; Moreau 
20081 1. This term is dominant for large Bo, and it depends on the second vertical deriva¬ 


tive of the velocity: for infinite Bq the bulk velocity field must be independent of the 
vertical. 

For finite but strong external magnetic field Bq, one may ask the following ques¬ 
tions: does the bulk flow remain three-dimensional, with weak three-dimensional velocity 
fluctuations superposed on a strong two-dimensional flow, or does it become purely two- 
dimensional? How do the Ohmic and viscous dissipation rates compare for large Bq7 Is 
there a scaling law relating the vertical and horizontal gradients of the velocity field? Or 
is there a critical value of Bq above which the flow becomes completely independent of 
the vertical? 

These central questions are strongly related to the scaling of the energy dissipation in 
MHD turbulence. Indeed, three dimensional hydrodynamic turbulence displays a dissi¬ 
pation anomaly: the 3D flow of a fluid with kinematic viscosity v forced on a scale £ with 
typical velocity U dissipates a finite power per unit mass e ~ U 3 /£ as the Reynolds num¬ 
ber Ut/v goes to infinity. By contrast, body-forced two-dimensional flows do not present 
such a singular zero-viscosity limit: in 2D the dissipated power vanishes for vanishing 
viscosity, £ and U being kept constant. Because the MHD situation with strong exter¬ 
nal Bq is an intermediate situation between 2D and 3D turbulence, an open question is 
whether it displays a dissipation anomaly. 

In the laboratory, the influence of an external magnetic held on the bulk turbu¬ 
lence of a liquid metal was studied in various geometries, ranging from channel flows 
to confined turbulence I Alemanv et al. 19791 : Gallet et al. 20091 : Klein fc Potherat 20101 : 
Pothe ra t fe Klein 112014 ). These studies indicated a clear trend toward two-dimensionalization 
as the external held increases. Experiments were also performed in shallow layers of liquid 
metal subject a strong uniform mag netic held to maximize the two-dimensionalization 
( Sommerialll98fil : lGallet et al. 11201 111 . The corresponding hows displayed good agreement 
with the phenomenology of 2D turbulence. 

Numerical simulations at unit magnetic Prandtl number have also displayed partial 
two-dimensionalization, with the formation of large-sca le two-dimensional structures co¬ 
existing with small-scale t hree-d imensional huctuations ( Alexakisll201lt iBigot fe Galtier I 
2011 1. Zikanov fo Thess ( 1998tl focused their numerical study on MHD turbulence in the 
quasi-static limit Rm —> 0, and they observed three regimes: for weak applied held, the 
how is three-dimensional. For intermediate Bo, its behavior is intermittent, with long 
phases of quasi-2D motion interrupted by 3D bursts. For large B 0 , the system settles 
into a purely 2D state, the amplitude of the ^-dependent part of the velocity held de¬ 
creasing exponentially in time. Although the forcing procedure is somewhat unrealistic 
(they rescale the total energy of the large-scale Fourier modes to a constant value at 
each time step), their simulations provide a clean example of MHD turbulence reaching 
a purely 2D stationary state provided the external held is strong enough. On the theoret¬ 
ical side, they consider a given two-dimensional large-scale cellular how, and they show 
that it is absolutely stable with respect to verticall y-dependent pertur bations for large 


enough external held. In a subsequent paper ( Thess fc Zikanov 2007t h they study the 


linear stability of several other steady two-dimensional hows as a function of the strength 
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of the external field. For finite Reynolds number, both the elliptical vortex and the vor¬ 
tex sheet are stable provided the applied field is strong enough. Although enlightening, 
these examples do not give an unambiguous answer to the central question asked above, 
because they do not consider the general setup of body-forced turbulence. 

Insightful analytical res u lts on the behavio r of such body-forced MHD flows have been 
obtained bv IPotherat fe Alboussiere 1 ( 2003h in the vanishing -Rm approximation: using 
rigorous methods from mathematical analysis, these authors provide upper bounds on 
the dimension of the attractor of MHD turbulent flows. For strong external magnetic 
field, the bound follows the same scaling as the one obtained (with the same method) for 
purely 2D flows. This indicates that the behavior of such low- Rm MHD flows somewhat 
“resembles” that of 2D flows, but it does not prove that the flow really is 2D: for instance, 
the bound holds even when the forcing is three-dimensional, a situation for which the 
flow cannot be purely two-dimensional. 

In this paper we use rigorous analysis and estimates to answer the central question 
of exact two-dimensionalization. We consider an MHD flow driven by a horizontal body 
force and subject to a vertical magnetic field, in a domain that is periodic in the horizon¬ 
tal and bounded vertically by stress-free surfaces. Such boundary conditions are chosen 
to remove the boundary layers inherent to experimental realizations, and to allow for a 
purely two-dimensional velocity field. We stress the fact that the results presented in this 
paper carry over to the 3D periodic cube, as discussed in the final section. In the quasi¬ 
static (Rm —> 0) approximation, we prove that the flow becomes exactly 2D asymptoti¬ 
cally in time, regardless of the initial condition, provided the interaction parameter N is 
larger than a critical value N c (Re). We call this property absolute two-dimensionalization. 
We then consider the full-MHD system and study the evolution of vertically-dependent 
perturbations of infinitesimal amplitude. We show that such infinitesimal perturbations 
decrease provided the interaction parameter is larger than some Re-dependent threshold 
value, and the magnetic Reynolds number is lower some other Re-dependent threshold. 
We thus extend the domain in which we prove that the two-dimensional flow is an at¬ 
tractor, but at the price of considering infinitesimal 3D perturbations only. We therefore 
call this property linear two-dimensionalization. 

In section [2] we describe the MHD system and recall some bounds on the viscous 
energy dissipation in hydrodynamic turbulence. Section [3] contains the proof of absolute 
two-dimensionalization for quasi-static MHD (Rm —> 0). Sections 0] and [5] deal with 
linear two-dimensionalization for the full MHD equations. In section [6] we summarize 
the different results obtained in this paper and discuss their implications for the energy 
dissipation rate of low-Rro MHD turbulence. Some technical aspects of the proofs are 
presented in three appendices. 

The rigorous estimates of this study are based on Young’s, Holder’s, Poincare’s, the 
Cauchy-Schwarz, and Ladyzhenskaya’s inequalities. Young’s inequality is simply the fact 
that \ab\ < -^a 2 + | b 2 for any c > 0, where c can be optimized to get the best possible 
upper bound. Holder’s, Poin care’s and the Cau chy-Schwarz inequalities are presented in 
a fluid dyn amics context in Doering fc Foiasl ( 2002 1. while Ladyzhenskaya’s inequality 
is proven in lLadvzhenskava' (jl963| i (Lemma 1. page 8). 


2. Body-forced magnetohydrodynamics and hydrodynamics 

2.1. Magnetohydrodynamics in a cubic domain 

We consider the setup sketched in figure [U an incompressible fluid of kinematic viscosity 
v, density p and electrical conductivity a flows inside a domain (x,y,z) £ V = [0 ,L\x 
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B x = 0 ,B y = 0 , d z B z = 0 



L 


Figure 1 . Flow of an electrically-conducting newtonian fluid subject to a uniform vertical 
magnetic field Bo. It is driven by a horizontal body-force f that is independent of the vertical. 
We assume periodic boundary-conditions in the horizontal, and stress-free boundaries with large 
magnetic-permeability at 0 = 0 and z = H. 


[0, L\ x [0,iT] with a Cartesian frame (e x , e y , e z ). It is stirred by a time-independent 
divergence-free two-dimensional horizontal body-force f(x, y) = (f x , f y , 0) that is periodic 
on a scale £, an integer fraction of L. That is, f(x,y) = |), where </> is periodic of 

period 1 in each dimensionless variable and has an rms magnitude 1. We refer to F as 
the amplitude, and </> as the shape of the force. 

A homogeneous steady vertical magnetic field £?o e z is applied to the system. We con¬ 
sider periodic boundary conditions in the horizontal directions (x and y), and stress-free 
boundary conditions for the velocity field at the top and bottom boundaries. We further 
limit attention to the case that these two boundaries have very large magnetic permeabil¬ 
ity, so that the magnetic field is vertical at z = 0 and z = H. Denoting u(x, y , z, t) and 
B (x, y, z, t ) as the velocity and magnetic fields inside the fluid, the boundary conditions 
at the top and bottom boundaries of the domain are 


0 , dgUy - 0 , U Z - 0 

B x = 0 , By = 0 , d z B z = 0 


| at z = 0 and z = H . 


( 2 . 1 ) 


The magnetic field inside the fluid follows the induction equation, which using the 
decomposition B(x, y, z, t) = Boe z + b(x, y, z, t) reads 


<9 f b + (u • V)b - (b • V)u = B 0 d z u + ? ? Ab , V.b = 0 (2.2) 


where rj = l/(/xocr), with fio the magnetic permeability of vacuum. Without loss of 
generality, we assume Bq > 0 in the following. 

The velocity field follows the body-forced Navier-Stokes equation with the Lorentz 
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force 

<9 t u+(u-V)u = -Vp+i/Au + f(a;,!/) + — (V x b) x (B 0 e z +h) , V.u = 0 . (2.3) 

PM o 

From equations (12.21) and (12.31) we define the kinetic and magnetic Reynolds numbers 
using the root-mean-square velocity U, where the mean is performed over space and time: 

Ul „ U£ . . 

Re = — , Rm = — . (2-4) 

v rj 


In the following, we focus on weak magnetic Reynolds numbers, Rm < 1, so that the 
dominant balance in (12.21) is 


R U t 

B ° £ P ' 


(2.5) 


One can define a third independent dimensionless quantity, the interaction parameter N, 
to evaluate the ratio between the Lorentz force and the inertial term in the Navier-Stokes 
equation: 


N = 


ctBqI 


Bli 


( 2 . 6 ) 


pU ppopU 

N is a dimensionless measure of the strength of the external field Bq. 

This magnetohydrodynamic problem admits two-dimensional solutions of the form 
u = V(x, y, t), b = 0, where V(x, y, t) satisfies a standard non-magnetic two-dimensional 
Navier-Stokes equation. Our goal is to prove that, for large enough interaction parameter, 
such motion is robust with respect to three-dimensional perturbations, i.e., the three- 
dimensional part of the velocity field decays in the long-time limit, and the system relaxes 
to a 2D flow. 


2.2. The quasi-static, or vanishing Rm limit 

The quasi-static limit consists in going to the limit Rm —> 0 while the interaction pa¬ 
rameter N remains of order unity, before expanding the equations to lowest order in 
Rm. At this order, the induction equation (12.21) reduces to pAh = —B 0 d z u, i.e. the 
induced magnetic field follows the velocity field adiabatically. It is much weaker than 
the external field, with b ~ RmBo- To lowest order in Rm, the Lorentz force is then 
(V x b) x B 0 e z /(pp 0 ), and the Navier-Stokes equation reads 

d t u + (u • V)u = — Vp + i/Au + f(x, y) — j3/±T 1 d zz u , V.u = 0, (2.7) 

where A -1 is the inverse Laplacian operator, which can be thought of as a division by 
k 2 of each Fourier component, k denoting the wavenumber. /? is the inverse Joule time 
scale corresponding to magnetic damping, /3 = crB^/p = Rg/(p/ioP)- It is related to 
the interaction parameter through N = [31/U. The main convenience of the quasi-static 
equation (12.71) is that it involves the velocity field only. 

2.3. The energy dissipation of 2D and 3D turbulence 
As mentioned in the introduction, the central challenge of the analysis is to establish the 
stability of some vertically-invariant flow \{x,y,t) with respect to vertically-dependent 
perturbations v(x, y, z, t ), without knowing the details of the base flow V: it is typically 
a turbulent solution of a two-dimensional body-forced Navier-Stokes system. Instead of a 
precise knowledge of V(a:, y, t ), for the proofs we make use of rigorous bounds on the time- 
averaged energy dissipation for solutions of the body-forced 2D and 3D Navier-Stokes 
equations. The corresponding bounds on the gradients of V are sufficient to establish its 
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stability with respect to 3D perturbations. We now recall these bounds, committing the 
proofs to an appendix. _ 


For hydrodynamic turbulence in a 3D periodic domain, iDoering fe Foias (1200211 de¬ 


rived a bound on the time-averaged viscous energy dissipation of any solution u to the 
body-forced Navier-Stokes equation, i.e., equation (12.31) without the Lorentz force. The 
method consists in bounding the amplitude of the body-force, and the power it inputs, in 
terms of the root-mean-square velocity U. We prove in appendix 1X1 that the bound easily 
carries over to the present magnetohydrodynamic system in the quasi-static, limit, i.e., 
when it is governed by (12771) . In the following, the constants with subscripts (cj)jgN are 
dimensionless positive numbers that depend on the shape of the forcing function only, not 
on its amplitude or length scale Z. By contrast, the undecorated c denotes the generic 
positive 0(1) constant that does not depend on any parameter of the problem. With 
these notations, the bound reads 


|Vu||!>< 


u 3 l 2 h 


(2L 


Zv \Re 


+ c 2 


( 2 . 8 ) 


where || ... || 2 is the standard L 2 norm over V, (• ■ •) denotes time-average, and Ci and c 2 
depend only on the shape of the forcing function f. 

Bounds can also be obtained for solutions of the 2D Navier-Stokes equation, i.e., for 
vertically-in variant horizontal vel o city fi elds satisfying equation (12.31) . Denoting the vor- 
ticity as u>, lAlexakis fe Doering ( 2006 1 derived a bound on time-averaged dissipation 
rate of enstrophy ||w || 2 in 2D. For a time-independent forcing function, and using our 
notations, their equation (19) translates into 


<l|Vcc||!)< 


hl 2 u 3 

vZ 3 


, c 3 


C4_\ 

Re) 


(2.9) 

( 2 . 10 ) 


where C 3 and C 4 are again dimensionless constants that depend only on the shape of f, 
and the prefactor H arises from || ... || 2 being the L 2 norm in 3D. 

This bound for vertically-invariant flows can be further reduced if the forcing is such 
that V 2 f is proportional to f. These harmonic forcings are sometimes called “single¬ 
mode” , or Kolmogorov forcings. For such forcings an improved bound on the enstrophy 
dissipation rate is 

u 2 l 2 h 

< 11 ^ 111 ) < c—^— . ( 2 . 11 ) 


3. Absolute two-dimensionalization for vanishing magnetic Reynolds 
number 

3.1. Evolution equation for the kinetic energy of the three-dimensional part of the 

velocity field 

We consider the evolution of an initially 3D velocity field under equation (12.71) . Define 
the vertical average of a field h as 

1 f z=H 

h (x,y,t) = — h(x, y, z, t) dz , (3.1) 

Jn J z -0 

and decompose the velocity field into 

u = V(a;, y, t ) + v(a;, y, z, t), with V = u, and v = u — V. 


(3.2) 
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In the following we call V (resp. v) the two-dimensional (resp. three-dimensional) part 
of the velocity field. Our goal is to find a criterion that insures the decay of v. One 
can obtain the evolution equation for the two-dimensional part of the velocity field by 
vertically-averaging the Navier-Stokes equation m, 

d t V + (V • V)V + (v- V)v = Vp + i/AV + f(x, y) . (3.3) 

By subtracting this equation to the Navier-Stokes equation, we get the evolution equation 
for the three-dimensional part of the flow: 

5 t v+(V-V)v+(v-V)V+(v-V)v— (v- V)v = -V(p-p) + v&v-0A~ 1 d zz v . (3.4) 

The incompressibility constraint imposes V • V = 0 and V • v = 0. To get the evolution 
equation for the kinetic energy of the 3D part of the velocity field, one takes the dot 
product of (|3.4p with v and integrates over the domain. Some integrations by parts using 
the incompressibility constraint lead to 

^ (3-5) 

with R a time-dependent quadratic functional of the field v, 

R{v,t} = - [ v(VV)-vd 3 x-^||Vv||2-/3||V- 1 5,v||2, (3.6) 

JV 

where we introduced the notation V -1 = VA^ 1 . 

The strategy is the following: we seek a function Ai it) such that 

• For all divergence-free fields v with zero vertical average (v = 0) satisfying the 
boundary conditions m, 

^{v,t} < Ar^llvlll. (3.7) 

• The long-time average of Ai is negative, i.e. (Ai) < 0. 

Then the kinetic energy in the three-dimensional part of the velocity field decays at 
least exponentially in time. The velocity field then becomes two-dimensional asymptot¬ 
ically in the long-time limit, regardless of the amplitude of the initial three-dimensional 
perturbation v. In the following, we call this property of the system absolute two- 
dimensionalization. Indeed, our approach for proving two-dimensionalization is very sim¬ 
ilar to proving absolute stability of some velocity field V to perturbations v. However, 
it is conceptually distinct from a standard absolute stability analysis: in relation (13.61) . 
V(a;, y, t) is a possibly turbulent velocity field that satisfies equation (13.31) . which includes 
an undetermined Reynolds stress term. 


3.2. A sufficient condition for two-dimensionalization 

To find a function Ai that satisfies (1X71) . we need to bound the first term on the right- 
hand-side of m - This task is performed in appendix [Bj where we use Ladyzhenskaya’s 
inequality on the L 4 norm. We obtain 




(VV) • vd 3 x 


/ IIVVII 2 

V lVh 


iiwiin 

vH ) 


IMH + ^||Vv|| 


2 

2 • 


(3.8) 


We bound the last term in (13.61) using Poincare’s inequality in the z direction only, 




—/?||d 2 V _1 v||| 


< -Pjp H V_lv ll 


2 

2 5 


(3.9) 
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to get the following bound on the functional "H, 


7£{v, t} < c 


piwil 


2 + II vv 


V ll 2 — / 3-^2 II V v ||2 — —|| Vv||, . 


\ LVH vH ' "" 2 

At this stage it is convenient to decompose v into Fourier components, 


v = ^v k (t)e d< 


+ c.c.. 


(3.10) 


(3.11) 


where the wavevector k takes the values 

( 27T 27T 7T 

k = (k x ,k y ,k z ) = [n x z j-,n y — ,n z — ) .with ( n x ,n y ,n z ) eNxNxM*. 


(3.12) 


The value n z = 0 is excluded because v is the three-dimensional part of the velocity field 
and therefore it has zero vertical average. 

Using this Fourier decomposition together with Parseval’s equality, equation (13.101) 
gives 




l v k| 


< 


/ ||VV || 2 [ IIwill 


V lVh 


vH 


|2X • ( P* 2 

-^Amp + 2 k 


(3.13) 


which we write as 

%{v,f} < Ai( 


1 2 , with Ai(£) = c 


p|VV||2 

V LVH 


IVVIII 

i/H 


2 ^- ( 3 - 14 ) 


To prove the decay of the three-dimensional part of the velocity field in the long-time 
limit, we need the time-average of Ai t o be negat ive. Using the concavity of the square- 
root function, we have (||VV|| 2 ) < y/(|| WHI), hence the following upper bound on 

<Ai>, 

<Ar> < c + ^ <11 Will)) - ^VW». (3.15) 

If the right-hand side expression is negative, then so is (Ai) and the flow becomes exactly 
2D in the long-time limit. This expression is a quadratic polynomial in \J (|| W|||). From 
the computation of the roots of the polynomial, we obtain that it is negative if 

<l|W|| 2 )<cy/^ 3 / 2 , (3.16) 


where c is yet another dimensionless positive constant that is independent of the param¬ 
eters of the problem. 

The term (|| VV|||) can be bounded in terms of the viscous energy dissipation per unit 
mass e inside the flow, 

<I|VV||^><<||Vu|© = ^<UA^(|L + c 2 ) , (3.17) 

where we used the bound (12.81) on the energy dissipation of a 3D flow. We obtain a 
sufficient criterion for two-dimensionalization expressed in terms of the root-mean-square 
velocity by demanding that the right-hand side of (13.161) be greater than the right-hand- 
side of (13.171) : 
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This criterion shows unambiguously that, whatever the Reynolds number is, there is 
a critical value of the interaction parameter above which the three-dimensional part of 
the velocity field decays and the flow becomes two-dimensional in the long-time limit. 
This two-dimensionalization occurs for three-dimensional perturbations v of arbitrary 
amplitude, provided this initial condition leads an rms velocity satisfying (13.18L and we 
therefore call it absolute two-dimensionalization. 

As a side note, we stress the difference between specifying the rms velocity U and 
specifying the amplitude of the forcing. All the initial conditions leading to an rms 
velocity U satisfying (13.1 81) two-dimensionalize in the long-time limit. By contrast, for 
constant amplitude of the forcing, two different initial conditions may end up in different 
attractors of the system, with differing rms velocities U\ and t/ 2 . It may be that U\ 
satisfies (13.1811 but U 2 does not. 


4. Linear two-dimensionalization 

In this section we come back to the full-MHD problem, with the goal of proving two- 
dimensionalization for finite values of Rm. We show that, starting from an infinitesimal 
three-dimensional velocity perturbation, the flow becomes two-dimensional asymptoti¬ 
cally in time provided the interaction parameter is above some threshold value and the 
magnetic Reynolds number is lower than some other threshold value. 

Let us decompose the velocity field into u = V(x, y , t ) + v(x, y, z, t), where V(x, y, t ) 
is a solution of the 2D body-forced Navier-Stokes equation, and assume that v(:r, y, z, t ) 
and b (x,y,z,t) are infinitesimal perturbations to this (possibly turbulent) base state. 
The three fields V, v and b are divergence-free. The linearized evolution equations for v 
and b are 


<9 4 v + (v • 

V)V + (V-V)v = - 

- Vp + zAW -I—— B 0 d z b , 

(4.1) 

PP 0 

d t b + (V 

• V)b = B 0 d z v + (b ■ 

■ V)V + yAh , 

(4.2) 


where p is the pressure perturbation. From these equations, we compute evolution equa¬ 
tions for the three following quadratic quantities 


d_ 

dt 


d_ 

df 

d_ 

dt 





- / v • (VV) • vd 3 x — Vv ||2 + —/ v-a,bd 3 x, 

Jv PPo Jv 

I b-(VV)-bd 3 x-?7 ||Vb||2-— / va z bd 3 x, 
Jv PP 0 Jv 


(4.3) 

(4.4) 


/ 9 z v-(VV)-bd 3 x- / b ■ (VV) • <9 z vd 3 x 
I'D JD 


+(z/ + 7?) f (Vd,v) • (Vb) d 3 x + B i 
JD 


(— 

Wo 


WdMl 


(4.5) 



With the goal of determining a Lyapunov functional for the system of linearized equa¬ 
tions, let us define 


•Fa{v,b} 



l|b||! 


+ a 


v • d z h d 3 x, 


2 2ppo 


v 


(4.6) 
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where a is a positive real number, and combine (14.311 . (14.411 and (14.5(1 to obtain 


d t 


P a = - v • (VV) • vd 3 x — 


|Vv||i + 


1 


b ■ (VV) • bd 3 x (4.7) 

b ■ (VV) • 9 z vd 3 x 

i 

+{y + rf) I (V<9 z v) • (Vb) d 3 x + B 0 ( — ||6»^b||| — ||<9 z v|| 

Jv \PH o 


J D 

— I 

PH 0 


Vbll 


PH o Jv 
d z v ■ (VV) ■ b d 3 x — 


IJV 


Because the base state is independent of the vertical and equations G3D and dH are 
linear, we can decompose the perturbation into Fourier modes in z and study the inde¬ 
pendent evolution of a single vertical Fourier mode. In the following, we therefore consider 
that (v, b) is a single vertical Fourier mode of perturbation, with vertical wavenumber q. 

Then ||<9 z v||| = O' 2 1|v||| and ||<9 z b||| = < 7 2 II b|||. We denote as Pa} the functional P a for 
such a Fourier mode, where the value of a q is free to depend on q. Using Holder’s, the 
Cauchy-Schwarz and Young’s inequalities, we obtain from m 


^ < IIWIIc 


l v ll 2 - HIVvl 


Ib||| 


llVVIloo^-lli--||Vb 


PH o 


V_ 

PHo 


(4.8) 


l|W|| 

1 + Pm 


\[Pm 


1 

n, 

H 

7T 


■a q q 2 \\h\\l 


+ B 0 q 2 -||b|| 2 - 

PHo 


y/m- H Vv ||^ + 


1 


T]7T 


2 pfioH 2 


|Vb|| 


where we used Young’s inequality (together with an optimization) to get the terms in 
the parentheses of the second and third lines, and we introduced the magnetic Prandtl 
number Pm = v/r] to alleviate somewhat the equations. 

We then decompose the Ohmic dissipation term into 


and rewrite (14.81) as 


— I 

PHo 


Vb| 


< — 
2 — 


2pfi 0 


IVbl 


2pho 


2 

2 ^ 


(4.9) 


d_ 
d t 


<Ml[2\\VV\\ x -a q B 0 q 2 


|Vv|| 


+ 


Nl r 


PHo 


IIVV||oo(l + a g q PHo) + a q B 0 q 


+ v- 


|Vb||g r 

PHo 


1 a q ir 2 ^pp^(l + Pm) 

2 2 H 2 VPm 


_ 1 a q g 2 H + Pm) 

2ny/ Pm 

2 _ m 2 ' 

2 


4.10) 


The method is now similar to the one used in section 13.11 if we can find a value of a q 
such that 

• is a positive functional, 

• there exists a real function \ 2 (t) such that < A 2 (t)P { a q J, 

• the time average of A 2 (t) is strictly negative: (A 2 (t)) < 0, 

then the vertical Fourier mode q of the perturbation decays to zero in the long-time limit. 
And if we can find such an ct q for all Fourier modes with nonzero vertical wave number q 
allowed in the system, then the L 2 norm of the ^-dependent perturbation decays to zero 
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in the long-time limit, and the system relaxes to the 2D solution of the Navier-Stokes 
equation. 

To find a suitable a q , we first impose the constraint 


7 r y/ Pm 

' _ 

9 q 2 Hjm{ 1 + Pm) ' 


(4.11) 


The reason for this constraint is made clear before equation (14.151) : it ensures that the 
terms proportional to || Vv||| and || Vb||| in (14.101) are negative and can be discarded. 

This constraint, together with q > it/H, Cauchy-Schwarz and Young’s inequalities, 
leads to 


QLq 


v • ebb d 3 x 


/x> 


< a q \ 


v 2 


(4.12) 


< a q \q\y/pjl^ ( 
\J Pm 


(\\nl 


2 ppo/ 


< 


1 + Pm 


+ 


2pp 0 J 


so that an upper and a lower bound on the functional Pa] are 


0 < 


2 pp, 0 J 


1 + Pm 1 ° 


l|b||| \ 

2ppo ) 


1 + 


V Pm 
1 + Pm 


or, since we do not care about optimizing the prefactors, 


»<i 


2 pp 0 




2 pp 0 J 


(4.13) 


(4.14) 


where we used the fact that -^p™ < Inequality (14.141) proves that PaJ is a positive 
functional. 

From the constraint (14.111) and using again q > ir/H, we show that the terms propor¬ 
tional to ||\7v11 2 and || Vb||| in (14.101) are negative, and upon dropping them we get 


-pM < 


dt 


a„ — II v II 2 


+ 1 


PM o L 
<4||VV|| 

Let us finally pick 


[2|| VV||oo — a q Boq 2 


II VV|| 00 (1 + a 2 q 2 pp 0 ) + a q B 0 q 2 - 


(4.15) 


1 - 


2pMo / 


a q B 0 q 2 \\v\\ 2 2 + a q B 0 q 


m 2 \ 

2 H 2 J 


PMo 


7 


OLn — 


B 0 q‘ 


(IIWIU). 


(4.16) 


(remember that we consider Bq > 0). We can do so if this value of a q satisfies the 
constraint (14.111) . Then 


df 


^ < 4||VV|| C 


l|b|||\ 

2ppo / 


^(IIVVIU) 


(4.17) 


+ 14(|| Wlloo)- 
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Let us further assume that we are in the region of parameter space where 


777T 2 

W 


>28(||VV|| 00 ), 


so that 


^4? < 4||VV|| 



^(livviu) (MI 



(4.18) 


(4.19) 


Using respectively the lower and upper bounds (14.141) on Fa} to bound the positive and 
negative terms in the right-hand-side of this inequality, we finally obtain 


^} < A 2 {t)F$> , with A 2 (t) = (silWIloo - y (IIVVIU)) (4.20) 


Because the function A 2 (t) has a negative time-average, the three-dimensional perturba¬ 
tion decays in the long-time limit. Of course, along the way to obtain (14.2011 we made 
some assumptions that constrain the region of parameter space in which the conditions 
hold. Indeed, we need the chosen value (14.161) of a q to satisfy the constraint (14.111) . and 
the magnetic diffusivity to be large enough for (14.181) to hold, i.e., we need respectively 


Bp 7(1 + Pm) (||VV||oo)g 2 tt 2 

(||VV||oo)Ry^ 7 r \f~Prn ,an V "28’ 


(4.21) 


They involve some kind of interaction parameter and magnetic Reynolds number built 
with the Loo norm of the velocity gradients. The former must be large enough and the 
latter small enough for magnetic damping to be efficient. Because one has little control 
over the infinite norm of the velocity gradients, we wish to express these constraints using 
the rms velocity only. 


5. The large Reynolds number regime 

We focus on the large -Re regime and therefore consider Re > 1. Then, for a solution 
V (x, y , t) of the body-forced two-dimensional Navier-Stokes equation, we can bound the 
time-averaged Loo norm of the velocity gradients (||W||oo) in terms of the enstrophy 
dissipation rate. This task is performed in appendix [Cl with the following result, 

(IIWIU) < (1 + lni?e + In , (5.1) 

where the H at the denominator arises from the fact that || ... H 2 is L 2 norm in 3D, and 
C 5 depends on the shape of the forcing only. The bounds (12.91) and (12.111) on (||Vw|||) 
then allow to bound (HWHoo) in terms of the rms velocity U. 


5.1. Time-independent forcing 

For a time-independent forcing, let us combine the upper bound on the vorticity gradients 
(12.91) with equation (15.11) to get 


(HWHoo) <C6 


1 + In Re + In — 


(5.2) 
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which we substitute into the inequalities ( 11 . 211 ) to get alternate sufficient criteria for 
two-dimensionalization, which now involve only the rms velocity of the base flow: 

2 


N > c 7 - 


H 2 L 2 


-Re 2 


and Rm < 




1 + In Re + In y 

cs 


H 2 L s/Re (l + In i?e + In j) ’ 


(5.3) 

(5.4) 


where we also assumed Re > 1 and Pm < 1 along the way, and C 7 and cs are positive 
numbers that depend on the forcing shape only. Forgetting about the logarithmic terms, 
the first condition requires the magnetic field to be strong enough for the interaction 
parameter to be larger than Re 2 . Note that the prefactor is proportional to H 2 , so that 
the constraint is more easily met in shallow domains. The second condition is to have 
a low enough magnetic Reynolds number, lower than roughly the inverse square-root of 
the Reynolds number. Once again, because the prefactor is inversely proportional to H 2 
this constraint is more easily satisfied in shallow domains. Condition (15.41) can be recast 
in terms of the kinetic Reynolds number using Rm = Pm Re , which gives 


Re < 


i 2 

(H 2 L)i 


_£s_ 

Pm (l + In Re + In j) 


2 

3 


(5.5) 


Neglecting the log term on the right-hand-side and assuming that the geometrical pref¬ 
actor is of order unity, the maximum Reynolds number for which we can prove two- 
dimensionalization is roughly Pm ~ 2 / 3 . For a liquid metal Pm ~ 10 -6 , so this is typically 
10 4 . Once again, in shallow domains (small H/£) two-dimensionalization can be proven 
up to higher Re. 


5.2. Kolmogorov forcing 

The tighter upper bound (12.111) on the time-averaged vorticity gradients holds if the 
forcing function is of single-mode type, that is if Af is proportional to f. For such forcing 
we combine (12.111) with equation (15.11) to get a bound on the velocity gradients in terms 
of the rms velocity 

(IIVVIU) < C 9 ^ (l + lnRe + lny) , (5.6) 

which we insert into (14.211) to obtain criteria for two-dimensionalization in terms of the 
rms velocity only: 


ip, / 

N > cio ——Re ( 1 + In Re + In — 




and Rm < 


c u 




H 2 L (l + lnf?e + In j) 


(5.7) 

(5.8) 


where we also assumed Re > 1 and Pm < 1 along the way, and Cio to Cn are numbers 
which depend on the shape of the forcing only. We can see that the assumption of single¬ 
mode forcing allows to prove two-dimensionalization in a larger region of parameter space. 
Indeed, forgetting about the logarithmic terms, conditions (EH) and (15.81) indicate that if 
the interaction parameter is larger than the Reynolds number and the magnetic Reynolds 
number is smaller than a constant, then the flow remains two-dimensional. Once again, 
the geometrical prefactors in (15.71) and (15.81) are such that the sufficient conditions for 
two-dimensionalization are more easily met in shallow domains. Condition (15.81) can be 
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System 

Region of parameter space 
(asymptotic in Re) 

Property 

Full MHD, time-independent forcing 

tt2 t 2 

N > const x — — — Re * 2 [ln(Re)] 2 
„ „ , l 3 * 

L2D 

f IL x bUIlo L A ,1 

H 2 L\PRe\n(Re) 

Full MHD, single-mode, time-independent forcing 

N > const x —-j—.Re [In (Re)] 2 
l 3 

Rm < const x „ . 

H 2 Lln(Re) 

L2D 

Quasi-static MHD, time-independent forcing 

N > const x ———Re 5 
£ 6 

A2D 


Table 1 . Regions of parameter space where we proved two-dimensionalization. The threshold 
values are written in the limit of large Reynolds number. L2D stands for “linear two-dimension- 
alization” and A2D for “absolute two-dimensonalization”. 


recast in terms of the kinetic Reynolds number only as 

Re<-^~ -4rr 7 -—-rr ■ (5.9) 

Pm H 2 L (l + lnRe + In y) 

For Kolmogorov forcing, we can therefore prove two-dimensionalization up to a Reynolds 

number of order Pm — . For a liquid metal, this is typically 10 6 . This maximum value of 
Re for which we can prove two-dimensionalization is larger in shallow domains because 
of the prefactor £ 3 /(H 2 L). 

6. Discussion 

In this study we considered a flow driven by a horizontal body-force and subject 
to a strong vertical magnetic field. We determined regions of parameter space where 
the two-dimensional solutions of the Navier-Stokes equation are stable with respect to 
three-dimensional perturbations. Such two-dimensionalization occurs for strong enough 
interaction parameter N and low enough magnetic Reynolds number Rm. We first fo¬ 
cused on the Rm <C 1 asymptotic limit, the so-called quasi-static limit, to prove absolute 
two-dimensionalization: above some Re-dependent value of the interaction parameter, 
the attractors of the flow are all 2D, i.e., the flow becomes 2D in the long-time limit, re¬ 
gardless of the initial condition. We then came back to the full-MHD problem, for which 
we derived a value N c (Re) of the interaction parameter and a value Rm c {Re) of the mag¬ 
netic Reynolds number such that, if N > N c and Rm < Rm c , the 2D flow is stable with 
respect to infinitesimal vertically-dependent perturbations: we called this property linear 
two-dimensionalization. Of course, the results proven using the full-MHD equations are 
also valid in the quasi-static limit Rm —> 0: then only N > N c is required for linear two- 
dimensionalization. Alternatively, one can directly prove linear two-dimensionalization 
for quasi-static MHD: bound the triple velocity product in (13.61) by || VVjlooHvH^, where 
V is now a solution of the 2D Navier-Stokes equation. The dissipative terms are bounded 
as in (13.131) . hence H{v, t} < (HVVjloo — y/2f3vTT/H)\\v\\\. Asking for the time-average 
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Figure 2. A sketch of the bounds (13.1811 . (15.71) and (15.91) . in parameter space, for single-mode 
forcing, i = L/ 2, and several values of the aspect ratio H/L. For simplicity, the dimensionless 
constants appearing in these expressions have been taken equal to unity. Top: Results obtained 
in the quasi-static limit. Above the linear two-dimensionalization line (blue, solid line), we have 
proven that the flow has a stable 2D attractor. However, this attractor is stable for small 3D 
perturbations only, and fully 3D attractors of the system may exist as well. Above the absolute 
two-dimensionalization line (red dashed), the flow has only 2D attractors. Bottom: Linear 
two-dimensionalization for the full-MHD system, for magnetic Prandtl number Pm = 10 -6 
(blue solid) and Pm = 10~ s (red dashed): on the top left part of each curve, we have proven 
that the system has a purely two-dimensional attractor. 
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of the parenthesis to be negative leads to the condition —y > for linear two- 

dimensionalization. Conditions expressed in terms of the rms velocity U are obtained 
using (ED followed by (12.91) or (12.111) : this gives the same criteria as (15.31) and (15.71) . 
albeit with different numerical constants. 


The different regions of two-dimensionalization are summarized in Tablc[T| and sketched 
in figure [2] in the case of single-mode forcing. These results have implications for the ex¬ 
istence of an energy dissipation anomaly in this system. Indeed, as can be seen in top 
panel of figure [2} depending on how one takes the double limit (Re oo ,N oo), the 
quasi-static system may or may not exhibit a dissipation anomaly. If we take this limit 
with N > Re 5 , then in the long-time limit the flow settles to a purely 2D attractor, 
with vanishing Ohmic dissipation, and a time-averaged viscous dissipation that follows 
a laminar scaling, proportional to Rer 1 : there is no dissipation anomaly (on a purely 
mathematical basis, our proof of absolute two-dimensionalization shows that the 3D part 
of the flow decays in the L 2 sense only, and one could argue that the 3D velocity gradi¬ 
ents can remain large. However, physical intuition suggests that the 3D velocity gradients 
decay when the L 2 norm of the 3D velocity decays). 

In the region N > Re (In Re) 2 , we know that the system has a 2D attractor that is 
stable with respect to infinitesimal perturbations. If the initial condition is close enough 
to being 2D, the system will end up in this purely 2D state, with no energy dissipation 
anomaly. However, this attractor is not necessarily the only one of the dynamics: if the 
initial condition is strongly 3D, the flow may end up in a fully-3D attractor, and this 3D 
state may exhibit a dissipation anomaly. 


The limiting values of the interaction parameter derived in the present study are 
rather conservative: the best criterion for two-dimensionalization we produced is roughly 
N > ReQxiRe) 2 , when several examples of 2D flows become stable as soon as N is larger 
than som e .Re-independent value (see for instance the study of the strained elliptical 
vortex by iThess & Zikanov I ( 2007t) h It could be that the computed threshold values of 
N are not sharp, or that the 2D turbulent flows we considered are much more unstable 
to 3D perturbations than the smooth ideal 2D flows considered by previous authors. 


We now wish to comment on the shape of the domain: we see that all the threshold 
values of N in Table Q] are proportional to the squared aspect ratio H 2 / L 2 of the domain, 
and the threshold Rm are inversely proportional to it: the shallower the domain, the larger 
the region of parameter space in which we can prove two-dimensionalization, as can be 
seen in figure [2] That is because shallow domains prohibit the long-vertical-wavelength 
perturbations that are poorly affected by magnetic damping. 


Further comments are needed regarding the boundary conditions considered in this 
study. We first stress the fact that the results of this study carry over to 3D periodic do¬ 
mains (albeit with different dimensionless numerical constants), which are more fashion¬ 
able among numericists. As far as laboratory experiments are concerned, the simple mag¬ 
netic boundary conditions in (12.11) can be relevant: they correspond to boundaries with a 
much higher magnetic permeability than the fluid. As an example, electromagnets often 
contain iron pieces to channel the magnetic field. If liquid metal is placed in the gap of 
such an electromagnet, between two iron pieces, the boundaries have very large magnetic 
permeability provided the ferromagnetic material is not saturated: the magnetic bound¬ 
ary conditions then approach those in (EH). Recently, such large-magnetic-permeability 
boundaries have attracted attention in the dynamo community, because in some geome¬ 


tries they l ower the threshold for dy namo action (IGissinger et al. I2008t iGallet et al. 


20121 20131 : Herault, fc Petrelisl 2014f) . However, these are not the most common bound¬ 


aries: many experiments and astrophysical situations correspond to non-magnetic and 
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insulating boundaries. These alternate boundary conditions do not affect the results ob¬ 
tained here in the quasi-static approximation, because equation (12.71) still applies. Indeed, 
both the high-magnetic-permeability and the insulating boundary conditions impose a 
vanishing normal component of the current density at the boundary, and therefore Neu¬ 
mann boundary conditions for the inverse Laplacian in equation (12.71) . As a result, both 
absolute and linear two-dimensionalization can be proven for insulating boundaries in 
the quasi-static limit, with the same criteria as for high-magnetic-permeability bound¬ 
aries. However, the proof of linear two-dimensionalization for the full-MHD equations 
does not carry over to non-magnetic boundary conditions, although we believe similar 
results probably hold. 

Finally, although the stress-free or 3D-periodic boundary conditions on the veloc¬ 
ity field might have some relevance in an astrophysical context, they are a crude ap¬ 
proximation to liquid metal experiments, where the fluid flows inside a rigid container. 
Such no-slip boundaries cannot be taken into account in the present analysis, because 


they prevent the existence of a purely 2D state. Instead, for strong enough magnetic 
field the flow settles in a quasi-2D state, with vertically invariant flow in the bulk 

and sharp Hartmann layers at the horizontal boundaries 

Sommeria & Moreau 1982: 

Klein & Potherat 2010: Potherat & Klein 

2014]). These layers can induce additional 

instabilities (Lingwood fc Alboussiere 1999; 

Moresco & Alboussiere 2004), and dissipa- 


tion through Hartmann friction, which we avoid by considering stress-free or periodic 
boundary conditions. 
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matics at UCLA. This research was also supported in part by NSF Award PHY-1205219 
(CRD) and by the junior grant TURBA from Labex PALM ANR-10-LABX-0039 (BG). 


Appendix A. Extension of the bound on time-averaged viscous 
energy dissipation to quasi-static MHD 

Doering fe Foias I ( 2002h produced an upper bound on the time-averaged viscous energy- 


dissipation rate for body-forced 3D turbulence in a periodic domain. We prove here that 
the same upper bound holds for MHD in the quasi-static limit. Following these authors, 
we write the forcing function as f = F<f>(x/£), where the shape-function <f> is periodic of 
period 1 in its two-dimensionless variables xjl and y/t. Its rms value is equal to unity. 
<f> is the shape of the force. 

To compute the kinetic energy budget, we take the dot product of (12.71) with u before 
integrating over the entire domain, which yields 


d (^ir) = X f ' ud3x _ Hl Vu|1 ^" • (A 1} 

After time-averaging, we obtain the following upper bound on the time-averaged viscous 
energy dissipation rate, 


v 



f- 



-PiWv-'dM 1 )<ful 2 h, 


(A 2) 


where we discarded the second term and used Cauchy-Scwhartz inequality to bound the 
first one. 
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To bound the forcing intensity F in terms of the rms velocity U, we take the dot prod¬ 
uct of (12.71) with f, integrate over the domain and time-average. The contribution from 
the Joule dissipative term vanishes after an integration by parts to move a ^-derivative 
onto the ^-independent force f. Further integrations by parts using the divergence-free 
constraint lead to 


|f|l!> = l 2 hf 2 < 


/ / u- (Vf)-ud 3 x) 

+ V 

/ / u • Afd 3 x\ 

\Jx> / 


\Jd / 


(A3) 


F 


< Cx-L'HIP* + c 2 - 


vFL 2 HU 


l 2 


where c.\ and £2 are dimensionless positive numbers that depend on the shape cj) of the 
forcing only. Multiplying by U/F yields 


r 2 rr 7 —it r - U 3 T , „ „vL 2 HU 2 

L~HFU < c\——L H + £2 -^- 


which wc use to bound the right-hand-si de of the inequality 
that is similar to Doering fc Foia s (1200 2). 


(A 4) 

. This gives a bound 


Il2 » _ U 3 r2rr „ vL 2 HU 2 
|Vu|| 2 ) < ci— L 2 H + c 2 


U 3 


i 2 


= L‘H~— I «. + K 


(A 5) 


Appendix B. Bounding the triple velocity product using 
Ladyzhenskaya’s inequality 

B.l. Ladyzhenskaya’s inequality in a 2D periodic box 

For this subsection, we consider functions of x and y only, and we denote as || ... || 2; n 
the standard L 2 norm in 2D. Ladyzhenskaya’s inequality provides a bound on the L 4 
norm of some field in terms of the L 2 norms of the field and its gradient. For a zero-mean 
scalar field f(x,y) in the periodic two-torus T 2 = [ 0 ,T] 2 , it reads 

JJ " r f 4 dxdy < C|| V/l||. J/||! ;X , (B 1) 

where C is a dimensionless constant of order unity. 

Let us consider a scalar field /(x, y) with a nonzero mean F over T 2 and a fluctuating 
part f(x,y), such that f(x,y) = F + f(x,y). We may bound the L 4 norm of f(x,y) 
according to 

[\ f^dxdy < F a L 4 + 8F 2 ||/|| 2 . ± + 3 [[ f 4 dxdy (B 2) 

JJr 2 ■> JTi 

< ^ll/ll2;J. +3 jj T f 4dxd y 

< ^ll/ll 2 ;n + 3C||V/|| 2 ;± ||/|| 2 ;± 

^ 3 (^ + ^) ll/l| A + ^ !||v/l| L. 


where we used Ladyzhenskaya’s inequality for the third step and Young’s inequality for 
the fourth step, with a some arbitrary positive number. Taking the square root of the 
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last inequality and using a convexity inequality, we obtain 

^ ff^ / 4 dxdy <y 3 + + L^W Vf\\% ± , (B 3) 

where the positive number a can be chosen to optimize the bound. 


B.2. Bound on the triple velocity product 

Let us write the three-dimensional part of the velocity field v = ( u,v,w ), and introduce 
Fourier series in z for the horizontal components u and v: 

p—-\- OO 

u(x,y,z,t)= ^2 u p{ x iV^) cos ’ (B4) 

P= 1 
p=-\-oo 

v(x,y,z,t)= ^2 v p (x,y,t) cos (pnjj) ■ (B5) 

P= 1 

(B6) 

The L 2 norm of one such horizontal component is ||u||| = y 2p=i~°° ll u pll 2 ; _L- 

Using these Fourier series in the vertical together with the Cauchy-Schwarz inequality, 
recalling that V is independent of z, we bound the triple velocity product according to 


ID 


TT P = + OO 

< Y I|W|| 2; x E 

P= 1 


V- (VV) • vd 3 x 

> 

and using inequality (IB 31) for u p and v pi we obtain 
v (VV) • vd 3 x 


JJ ujdxdy + J j J v^dxdy 


id 


< llVVIh: 


3 C 


3+^(IHl2;X+IMl2 !A ) 

+L^(\\Vu\\ 2 2 . x + \\Vv\\ 2 2 . ± ) 


< 


liwib 

Vh 


3 + 7^-|| v ll2 + L 


3Ca, 


Vvll 


,(B7) 


(B 8) 


In this expression the positive number a can be chosen arbitrarily to optimize the bound, 

2 TT # # # 

and we set a = ecXTjjVvlp' Using the concavity of the square-root function we finally 
obtain 


v (VV) • vd 3 x < 


ID 



Vv 


, (B 9) 


Appendix C. Bound on the time-average of the L°° norm of the 
velocity gradients 

Throughout this appendix we consider 2D functions of x and y only, and for brevity we 
denote ||... || 2 as the L 2 norm in 2D. We wish to bound the time average of the L°° norm 
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of the velocity gradients of a solution of the body-forced 2D Navier-Stokes equation. This 
bound should depend on the rms velocity only. 

Decompose V into Fourier components, 


V = ^V k (f)e lkx + c.c., (Cl) 

k 

where the wavevector k takes the values 

/ 27r 27r\ 

k = ( k x ,k y ) = \n x — 7 n y — j ,with ( n x ,n v ) <E N x N \ (0,0). (C2) 

Consider an arbitrary cut-off wavenumber Q > 0 and write 

IIWIU < £fc|v k | = E fc l v kl + E fc l y kl ( C3 ) 

k k<Q k>Q 

< +Jn 

y k<Q y k<Q y k>Q y k>Q 

Now the sums that are independent of V can be bounded by integrals, with the following 
result 


EpSi 2 ln<5, 


k<Q 


E 

k>Q 



fc 4 ~ Q 2 


(C 4) 


where the notation a < b means that there exists c > 0 such that a < cb 7 where c is 
a dimensionless positive constant that is independent of the parameters of the problem. 
Hence 

IIWIU^^HAwlla + v/i^QIIVwlla. (C5) 

Upon picking the value Q = L\\ Aw^/v^D, where D = IV w ll! + (II Voj|||) is a strictly 
positive time-dependent quantity - this value of D allows for the following simple upper 
bound on D~ x , D~ x < (||Vw|||) , which is needed to obtain inequality (1C 161) - and 

time-averaging, we obtain 

(IIWIU) < (Vd) + ^||V w || 2 y/h, (C6) 

and using Young’s and Jensen’s inequalities, 

(IIVVIU) < A /(||Va;|||) (l+hi^ 1,211 ^ 11 ' )) ■ (C7) 

We now bound the argument of the logarithm. Our goal is only to prove that it has at 
most an algebraic dependence on Reynolds number, so that the logarithmic term admits 
a bound of the form const x ln(i?e). The 2D Navier-Stokes equation written for the 
vorticity is 

a t u + V-Vw = I/A u + (V x f) • e 3 (C 8 ) 

Upon multiplying by - Aw and integrating over space, we get the palinstrophy evolution 
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equation 


d /||Vw|| 2 \ = /■ ( V . Vw)Awdx —i/||Aw|||- f (V x we z ) • Afdx (C9) 

di \ 2 / ■'[o.L] 2 ■/[0,L] 2 


Division by D yields 
Id 

2d t u J[ 0i l]2 

and after time-averaging 


(In D) = ^-l (V • Vu)Aw dx- 
D J [o.Ll 2 


.IIAwlli 1 


D D 


[ (V x we z ) • Afdx, (CIO) 
J[0,L] 2 


»-a. 


[0, Zy] 2 


(V- Vw)Awdx 


l|Aa;||| 

D 




F [0 , L ] 2 


(V x uie z ) • Afdx 


)' 

(Cll) 


To bound the first term, let us first use the Holder-Cauchy-Schwarz-Young recipe 


/ (V-Vw)Awdx 

/ [0 , L ] 2 

and divide by D to get 


<||V|| 00 ||V W || 2 ||A w || 2 <^||Aa;||2 + i- 


IVwlH, (C12) 


1 

D 



(V-Vu)Aw dx 


2D 



(C13) 


Since we only want to prove algebraic dependence on Reynolds number, the following 
simple estimate is sufficient at this stage: 

\\n 2 co S l|V|| 2 ||- v w || 2 < ^77^(||Vc|||)||V||I + "V ||Vo;||l. (C 14) 
ZLU 2 v(l|vw||^) 

The first inequality is a standard Agmon’s inequality for vector fields in a 2D periodic 
domain. After time averaging, we obtain (||V||;^) < UL\J (|| Vw|||). The time-average of 
(1C 131) gives 


1 

D 



(V- V«)Awdx 



(C15) 


Using Cauchy-Schwarz inequality, together with the convexity inequality (||Vu;|| 2 ) < 
>/I Vw|||) and Poincare’s inequality, the term involving the forcing is bounded by 


( -pr f (V x we z ) • Afdx 

\ V J [ 0 , L ] 2 


< 


D 


(C16) 


IIV^lbllAf|| 2 \ < IIAf|| 2 < L||Af|| 2 

(IIV^III) / - i)- mu 


Recall that we write the forcing function as f = F<f>{x./i), where <f> is periodic of period 
1 in its two dimensionless variables and has an rms value 1. Then following Alexakis and 
Doering, we can bound F in terms of the rms velocity U : their equation (18) gives 


F < 


IP 

T 



hence ||Af|| 2 < ^ (cs + 


(C17) 


where the (cj)igN appearing throughout this appendix are dimensionless positive numbers 
that depend on the shape 0 of the forcing only. 






























22 


B. Gallet and C. R. Doering 


Upon inserting (|C 151) . (1C 161) and (1C 171) into (1C 111) . we obtain 


L 2 ||Ac|| 
D 


< Re 


L 3 V( l|Vo;|||) , L 4 C/ 


Iv 


Pv 


C3 


Cj 

Re 




(C18) 


(C19) 


which is inserted into SC 71) to get an upper bound on the infinite norm of the velocity 
gradients in terms of the L 2 norm of the vorticity gradients 


(livviu^divciii) 


1 + In I 1 + Re 


£V<II^||^> 


£v 


+ ^4 (c 3 Re + c 4 ) 


(C 20) 

Since bounds on (|| Vw|||) ca n be computed in terms of t he rms velocity for body-forced 
two-dimensional turbulence ( Alexakis fe Doering 1 12006 1. the right-hand-side of (1C 201) 
can be expressed in terms of the rms velocity only. These bounds show that the argument 
of the logarithm has at most an algebraic dependence in Re, hence, for Re > 1, 


(C21) 

where £5 depends on the shape of the forcing only. This constant is denoted as C5 in the 
main body of this article, where the 2D L 2 norm || Vcj ||2 of this appendix is replaced by 
its 3D counterpart. 


(IIVVIU) < C 5 yj (II Vw||!) ( 1 + In Re + In j 
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